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Homogenization of a Boundary Obstacle Problem 

Ray Yang 



Abstract 



We prove the existence of a homogenization hmit for solutions of appro- 
priately formulated sequences of boundary obstacle problems for the Lapla- 
V^ , cian on C^'" domains. Specifically, we prove that the energy minimizers u^ of 

/ iVuepdx, subject to u > cp on a subset Se, converges weakly in H^ to a limit 
PLh ' ^ which minimizes the energy J \Vu\'^dx + Jj.{u — ip)'^fj,{x)dSx, S C dD, if the 

•^ , obstacle set /S^ shrinks in an appropriate way with the scaling parameter e. 

(-H I This is an extension of a result by Caffarelli and Mellet [1] , which in turn was 

"^ ■ an extension of a result of Cioranescu and Murat [3] . 

S: 

1 Introduction 

— ( ■ 

> 

Tij- I In this paper we consider the homogenization of a sequence of boundary obstacle 

problems. The boundary obstacle problem was used by Duvaut and Lions |1] to 
describe a semipermeable membrane, that is to say, a membrane which permits fluid 

i/^ . to flow in one direction, but not the other. We can think of the obstacle function 

O I ip[x) as representing the outside pressure or concentration, and S^ as the subset of 

the boundary which is composed of a semipermeable membrane that permits fluid to 
flow, Se being composed of "patches" of boundary, of appropriate size and separation 
given by e. 

It is natural to ask, as the scaling parameter e — )■ 0, what the limiting behavior 

5^ i might be. This is the realm of homogenization theory. Caffarelli and Mellet |T] 

answered this question in the case where S^ lies along a flat portion of the boundary, 
and consists of patches contained in small balls about points of a lattice eZ", with 
the capacity of each patch given by a stationary ergodic process. They found that, 
the minimizers m^ to the energy functionals 

J{u) = / |V«|^(ia; 
Jd 

converged to the energy minimizer of a special energy functional 

Ja{u) = / \Vu\'^dx + / a{u — LpY_dSx 
Jd Jt, 

1 



which is the regular energy plus an extra penalty term along S, the flat portion of the 
boundary. Caffarelli and Mellet's result derived from their earlier work on homoge- 
nization of an obstacle in the interior of the domain |2], and used the framework laid 
down by Cioranescu and Murat [3j. Focardi has analyzed these problems in the con- 
text of F-convergence and in terms of minimizing a nonlocal fractional order energy 
on the boundary surface [6] [5]. 

In this paper, we consider give conditions on the obstacle sets S^ that are sufficient 
for there to be a homogenization limit, dispensing with the lattice and allowing them 
to live on a portion of the boundary that is locally a C^'" graph. Specifically, we 
show that so long as a particular measure associated with Se converges to a surface 
measure in the appropriate space, the corresponding minimizers converge weakly to 
the minimizer of an energy with an added term. It is no longer required that each 
patch be centered on a lattice point (indeed, defining an appropriate lattice on a 
curved boundary is difficult), merely that each patch be well separated (a distance of 
about e) from every other patch, and be contained inside a small ball (specifically, of 

n-l 

radius Me 1-2 for some M). 

Our method of proof tracks closely with the framework followed by Caffarelli and 
Mellet [1]. In Section 2, we show that the homogenization result depends on the 
construction of an appropriate corrector, which is a family of auxiliary functions that 
captures the effect of the obstacle set on the solution. In section 3, we construct this 
corrector for the special case of an obstacle T^ consisting of balls centered on points 
of the boundary, and in section 4 we demonstrate that this corrector can be modified 
to the more general case of obstacles living on the boundary proper. 

2 The Problem 

2.1 Statement of the problem 

Consider an open bounded C^'" domain D C M", with n > 3, whose boundary is 
divided into two parts: 

aD = F U S 

where F fl S = 0. The first part, we call F. On the second portion, which we call S, 
there is a subset, known as the obstacle set, S^ C E. We consider the minimizers of 
the Dirichlet energy 



J(v) = / \Vv\'^dx 
Jd 



among functions v{x) satisfying Dirichlet conditions on F and a boundary obstacle 
problem on S^', we call this set the collection of admissible solutions to the obstacle 
problem: 

1. ve H\D) 



2. v\r = ip{x) for some smooth boundary data ip G H^{D) 
3- v\s^ > (/? for a smooth obstacle function ip G H^{D) 
There is a natural corresponding set of test functions, which satisfy 

1. (peC^iD) 

2. 0|r = O 

The existence of a solution to the boundary obstacle problem is assured by the 
usual considerations of Hilbert space theory. 

To prove our main theorem, we will rely on a model obstacle set, which does not 
live on the boundary, but close to it. We will then be able to extend our proof from 
this model set to the general case stated in our problem. This set is constructed as 
follows: 



Ts= fu^'-..(^^-'^))n^ 



where the x^^k are chosen to lie on S at a distance of at least 2e from each other, and 
the radii of the balls are chosen so that there exist constants Ci,C2, independent of 

k, e, so that 



n-l 



where Ci < f^ ^ < C2. Clearly, the number of such balls is 0(e^^"). To each obsta- 
cle set, we can assign a corresponding density function fieix) living on D, which is 
associated with 

k 
The model obstacle set can be generalized to arbitrary patches living on S, pro- 
vided certain conditions hold. It is sufficient that each patch be a set of appropriately 
bounded capacity, which is contained within B n-i (x^^k) H S, where the x^^k are 
again points distributed on S, at a distance of at least 2e from each other. We can 
call this set S'^^^ and its capacity is 7e,A;£"~^. The overall obstacle is then 



Se = \JSs,kC^CdD 



k 

The density corresponding to this is 

k 

This definition in the case of balls is in accord with the one previously given, up to a 
constant which depends only on the spatial dimension. 
We now state our main theorem: 



Theorem 1. So long as fi^dx — > fi{x)dSx in H^^ as e ^ 0, where fi{x) is some 
density on S, then the minimizers u^ -^ uq in H^{D) where uq is the minimizer of 
the energy 

^m(^) = / \^v\'^dx + Cn / {u - ipf_^{x)dSx 

J D Jt. 

among all functions in H^{D) satisfying the Dirichlet conditions on V. The constant 
Cn depends only on the spatial dimension n. 

Here we abuse notation slightly to take H^^ to be the dual space of iJ, defined as 
the subspace of H^{D) whose trace reduces to on T. 

2.2 Reduce Problem to finding a good corrector 

Following Caffarelli and Mellet, we introduce the notion of a corrector function, and 
demonstrate that the existence of a corrector meeting the appropriate conditions is 
sufficient to prove Theorem [H The proofs of this subsection are essentially unchanged 
from [Ij. 

Lemma 2. Assume that the obstacle set S^ on the boundary is as defined. Then for 
each e there exists a function We(x) G H^{D) such that the following conditions are 
satisfied 

We{x) = 1 if X e Se (1) 

||'y^'£||L°°(D) ^ C where C is bounded independent of e (2) 

u)e ^ zn H^{D) (3) 

and for every sequence Ve{x) satisfying 

f £ > for X & Se 

\\v£\\l°^{d) < C 

Vs -^ V in H^{D) 

and for any test function (f) G C°°{D) whose support is bounded away from T we have 

lim / Vws ■ Vve(f)dx > — / vcl)fi{x)dSx (4) 

^-^0 Jd Js 

with equality if Ve = on S^ . 

The functions {w^} are called the corrector, and demonstrating the existence of a 
corrector satisfying those conditions is the subject of most of the paper. Two auxiliary 
lemmas are needed before we prove that the existence of such a corrector is sufficient 
to prove our theorem. 



Lemma 3. For any test function </>€{(/)€ C°°{D) : 0|r = 0} we have 



lim / |Vwe| (j)dx = / (j)^{x)dS: 



^^0 'D JT. 



Proof of Lemma\^ Let v^ = 1 — w^ where w^ is the corrector. Then ti^ = on T^, 
f £ ^ 1 weakly in H^{D), and thus the resuh follows trivially from Lemma [2J D 

Lemma 4. Let u^ he a (suh) sequence of minimizers in H^{D). If Us ^ u weakly in 
H\D), then 

=^0 Jd Jd Jy. 



Proof of Lemma\^ For an arbitrary function v G H^{D) satisfying the Dirichlet con- 
ditions on r, we consider the function v + {v — ip)^Ws, which is admissible for the 
obstacle problem at scale e. We evaluate the quantity 

\Vu, - V(w + {v- ip)^We)\'^dx > 

D 

If we take the lim inf as e — > 0, we can take advantage of the fact that w^ ^^ in 
H^{D), whence we can write 



liminf / \Vue\ + \Vvf + \Vwe\ {v - v?)! 
Jd 

-2Vue ■ Vv - 2{v - (p)-Vu^ ■ Vwedx > 
We can apply Lemma [3] on the term in |Vi«£p to get 

/ |Vt«£p(f — ipf_dx — !■ / (v — LpY_^{x)dSx 

The term in V^e ■ Vwe can be rewritten through inequality (jl]), giving us 

liminf / Vue ■ Vwe{v — (p)-dx > — / {u — (p){v — (p)^fi{x)dSx- 

We are justified in invoking the lemma as u^ G L°°{D) (owing to the maximum 
principle and the boundedness of both the boundary data i/j and the obstacle function 
if), and because as the solution to an obstacle problem Us — ip > on the obstacle 
set Ss- 



We thus have 

hminf / Wu^l'^dx > - / \Vv\^ + 2Vu -Vvdx 

{{v - (pY_ + 2{u - ip){v - (/?)_) ^{x)dS^ 



D J D 



If we let t; — 7- -u strongly in H^{D) (and hence strongly in L^), then we can use the 
fact that {u — Lp){u — (p)^ = —{u — Lp)"^, getting the final result that 

'^^O Jjj J^ J^ 

D 

We now prove that the existence of a corrector is sufficient to demonstrate our 
theorem. 

Proof of TheoremUl Consider for the arbitrary test function v, the function v + {v — 
(f)-We, which is admissible as a possible solution for the obstacle problem at scale e. 
Applying Lemma [3] and expanding, we find that 

\imj{v + {v - (p)-We) = J^{v) 
However, we also know that 

since u^ is the energy minimizer of all admissible functions at scale e. Thus we have 

limsup J'(iie) < J'^(v) 
for any test function v. However, Lemma H] tells us that 

liminf J^(n,) >J^{u) 
where u can be any subsequential limit of the w^. Thus, 

JA^) < Mv) 

for arbitrary test functions v, which is to say any subsequential limit of the -u^ is the 
energy minimizer of J7^. Thus, the u^ converge to the energy minimizer of J7^. D 



3 Construction of a corrector 

We will first construct our corrector for the case of the model obstacle set T^ consists 

n-l ^ 

of arbitrarily placed balls of radius r^ ^ = r^^ke^^-^ centered on points x^^k ^ S? where 
the Xgr^k are have distance at least 2e from each other. The corrector for S^ will be an 
adaption of the corrector for T^. The construction of the corrector borrows directly 
from [3] , and it is nothing more than a truncation of the suitably rescaled fundamental 
solution for the Laplacian. Let 



We,k{,x) = 1 for |x - Xe,k\ < re,k 
1 / 1 






for r.fc < \x- Xe,k\ < £ 



•e.k 

for £ < \x — x^ k\ 



Each of the We,k is supported on a ball of radius e about some point Xs,k, and by 
assumption these balls are disjoint, so we simply take 



k 



as our corrector. 



3.1 Simple properties of the corrector 

That the corrector as constructed satisfies ([T]) and ([2]) is obvious. To demonstrate 
weak convergence in H^{D), or the condition in (^, we show that ||ti'£||L2 — )■ 0, 
and that ||Vwe||L2 < C where C is independent of e. These are two reasonably 
straightforward calculations. 



Lemma 5. 



£—>-0 ^ ' 



Proof of Lemma\^ We first perform the calculation inside a single ball of radius e. 
Without loss of generality, we assume that this ball is centered around 0, and we 



n— I 

recall that r^ = 0{e"-^). 



wJ'^dx = \BrA + I IwJ'^dx 



B^ Jb^\B. 



(4-n)(n-l) 2("-l) i o „. ^, "("-!) i .1 o„ . 

n(n— 1) on/ /t 

+ 0(£ "-2 ) + 0{e "-2 +^ ") + 0{e "-2 +^ ^"^ 

Since there are on the order of e^~^ such balls, our result holds. Note that we proved 
the result assuming the corrector is defined in the full ball. This is, of course, not the 
case, as the corrector is only defined in the intersection of D with the ball. However, 
the bound on the full ball is a good upper bound. D 



Lemma 6. 



|V«;^||^2(B) < C 



Proof of Lemma Again, we perform the calculation inside a single ball, as we did 
for Lemma [5l 

r ( |2;|i-" \^ 

iVwJ^dx = / \{2-n)—^^ — ,— \ dx 

' ^ ' / , \ ^ ' »,2— n p-2—n j 



2 />£ 

na(n) ( -r, ?; — / |a;|"'"^"(i|x| 

nain) ^f^^^" - r'^'"') 

^ ' ' ^2-n c-2-n^2 ^ "^ ' 



2-n 

2-n 

"-£ 
2-n 



na(n)^: — = Oie 



n-l\ 



c-2— n j,2—n 



Thus, since there are 0{e^ "") such balls, we have that the gradients of We are uni- 
formly bounded in L^. D 



3.2 Limiting property of the corrector 

Finally, to demonstrate (jl]) we rely strongly on fact that dD is locally the graph of a 
C^'" function. In what follows, we assume that e is sufficiently small. 

/ (j)VviNwedx = 2^ / (pVveVwedx 

SO we can examine the support of We about each of the Xs,k individually. Integrating 
by parts, we have that 



(pVv^VWedx = / (j)VedyWedSx 

(Be\Br^)r\D J d{(B,\Br^)r\D) 

since Aty^ = in the annulus. We separate the three portions of the boundary: 



+ / + 

d{{Be\Br^)UD) JdBr^nD J dDn{Be\Br^) J dBeHD 

The first integral is purely positive, so it can be neglected for purposes of proving 
the inequality. It is easy to see that it is when f^ = on T^, so it can also be 
neglected for proving the equality in that case. We will show the second integral goes 
to 0, and the third integral retrieves the desired term along the boundary in the limit. 

3.2.1 The Second Term 

Lemma 7. The effect of the second term, 

y / (pv^dyW^dSx — )■ 

as £ — )■ 0. 

Proof. We examine the second integral in a single ball. Here, we consider 

(j)VeduWedSx 
dDniB,\Br,) 

The idea is that dD is "almost fiat" for sufficiently small e, and that Ws being a 
purely radial function, its derivative in a transverse direction is 0, so its derivative in 
an "almost transverse" direction is small. 

Formally, consider a system of coordinates about x^^k, which for ease of notation 
we take to be 0, and where e„ represents the normal vector to the surface dD at 
0. The surface dD can be represented as x„ = S(x'), where x' are the first n — 1 
coordinates. Then E(x') is a C^'" function with tangent plane at 0; so specifically 
we have the estimate 

|S(x')| <C|xV+". 

9 



Let u{x') represent the normal vector at the point {x',Il{x'). This is a C" function. 
Now we can recast the integral as 



Now we can write 



f (pv.Vweix', J:{x')) ■ iy{x') ^/l+W^i^dx' 

J Br\Brcr 



(5) 



Vwe{x', S(x')) ■ iy{x') = Vwe{x', 0) ■ //(O) 

+ {Vweix', i:{x')) - Vw,{x', 0)) • Z/(0) 
+Vw,{x' ,T{x')) ■ {u{x') - u{0)) 

The leading term here is by construction. We will bound the second term by the 
(ji,a pYopeitj of S, and the third term with the C° property of u. We first obtain a 
pointwise bound on the components of Vwe{x', 0) — Vwe{x' , S(a:')) which are not in 
the Xry direction: 



.^^^(^',S(x'))-?^fa;' 



dxi 



dxj 



< 



0<a'„<S(x'') OXiOX. 



< — 



1 



sup 



{x',Xn)mx')\ 

n(n -2)\xiXn\ 



< C„"^""'^ .'^;jx' 



/|2+2o 



Sfx') 



j,2—n p2— n |^/|7i+2 ' 



< c 



nin — 2) 



j,2—n c-2— n 



\X 



l\l+2ct-n 



The bound on the derivative in the x„ direction is similar, but includes an extra term: 

^2„ 



|^(x',S(x'))-^(x',0)| < sup |^(x',xJ||S(x') 

dXn OXn 



Q<x„<T.(x') dx^ 



< 



sup 



n(n-2)xl 2- 



\n+2 



+ 



n 



< C ."~V ( 



j,2—n p2—n 



n\x 



/|l+3«— n 1^ I /ll+o— n 



\Ux') 



+ x' 



) 



We recall that the regularity of the boundary lets us pick e sufficiently small so that 
|VS| < 1 for |x'| < £, and also that 

\u{x')-iy{0)\ < C|x'|". 

Thus, equation \5\ can be directly bounded: 



10 



2\\(f)v,\Lo^ f Vweix', 0) ■ 1^(0) + (Vweix, S(x')) - Vw,{x\ 0)) ■ u{0) 
+Vw,{x', J:{x')) ■ {u{x') - u{0))dx' 

< 2\\(PVe\\L^ f \VWeix', S(x )) - VWe{x' , 0) | + \VWe{x' , E{x'))\{C\x'\'')dx' 



C 


-"J 


r2-« - £2- 

c 


-"J 


r|-" - £2- 

p.n+«— 1 


-n 



're 

c 



< C\\(l)V,llLo-^_^^ 



Since there are on the order of e'^ " such balls k on the surface, the total integral 
from this term is 0{e°'), which goes to as e —t- 0. D 

3.2.2 The third term 

With only one term remaining, we see that to prove Lemma [2] it suffices to show 
Lemma 8. As e — t- 0, we have 

lim y^ / (j)VeduWedSx ^ ~ / (j)viJ,{x)dSx (6) 

with equality if v^ = on the T^. 

Proof. As before, we start by examining the behavior confined to a single ball, which 
for ease of notation we assume to be centered at 0. Notice that along dB^, we have 



^ly^e ^T^e ^2-n _ ^2-n\-^\ 'e X _ ^^n-2 " 

which is thus constant over the surface of the ball. 

This inspires us ( following Cioranescu and Murat [3] ) to introduce the auxiliary 
function 

,x i^:^^^^H\x?-e^) for |a;|<£ 

q^[x) = Vl "sj 2eVi I ; II- 

for Ixl > £ 



11 



Thus, 

A 1 — n n , , 



r^" - e e 



while the normal derivative along the boundary of dBf. matches that of the corrector 

We- 

We see that ||q'£||L°o < Ce, and similarly that 

l|Vg.|li.(^^) < Ce- 

for every 1 < p < oo. If we repeat the construction of q^ over all 0{£^^'"') such balls, 
then we have 

||Vg,||^,(^) < Ce 

which goes to as e — > 0. In particular, q^^ — )■ in H^{D). Integrating over all balls, 
we have 

Jd Jd ^ Jd{B,{xk)nD) 

The second term on the left hand side is of order e and so can be neglected, while the 
term on the right is the sum of the term in ([6]), and another error: along dD fl B^, 
we have Vg^ = .^Jj^ -x. However, dD is a C^'" surface - using local coordinates 
and letting the tangent plane at be {x„ = 0}, we see that \xn\ < Clx'l^"*"", while 
the normal vector to dD is a C" function, with \u{x') — i^{0)\ < C|x'|", where x' 
represents the other n — 1 coordinates to x, and z/(0) = af„ by definition. Thus 

,^ , , 1-n 1 

Vge ■ z^ = \— x ■ ul 

r^^" — e e 

Integrating this over a surface of area 0{e"'~^), where we have 0(£^~") such surfaces, 
gives us a contribution of 0{e°'), which goes to as £ vanishes; in other words, the 
integral along dD is, up to a small error, negligible, since dD does not differ much 
from the tangent plane to which Vg^ is orthogonal. 

Hence, the lemma simplifies to evaluating the effect of 



Z^ ^2-n 



n I I IT- , 

(pVs—dx 






But we see that ^i},^_ = r"^^^ 1„_^ , so it can be bounded from above by r^ f, j^t^, 

and from below by r"^^. Thus we see that we have reduced ourselves to evaluating 
the effect of 

(f)VefJ.e{x)dx 
D 

12 



which by hypothesis converges to 

— / (f)viJ,{x)dSx 
as desired. D 



4 Adapting the corrector to randomly shaped ob- 
stacles 

4.1 The adapted corrector 

In this section, our goal will be to demonstrate that we can still fulfill the corrector 
conditions given in §2.21 for obstacle sets that are not composed wholly of small balls. 
To be precise, we replace Br^f,{xk) with an arbitrary set Se,k C S n-i {xe,k) H E, so 
long as this set has a well-defined capacity equal to 7(6, A;)£"~^, where 7 is bounded 
from above. The capacity is defined with respect to R"'. 



7(£:, /i;)e" ^= sup / \Vv\^dx 



The method of this section, which we borrow from [2], is to build a corrector by 
"stitching" together a translation of the capacitary potential of Ss k inside B n-i (xk) 

' Me"--' 

with the corrector we have already constructed for some "equivalent" radius r^ ^ up 
to Be{Xk). 

Define r^^k to be the radius of the ball with capacity 7(5, k)e"'~^. Since the capacity 
of a ball of radius r in R" is Cnr^'"^ where c„ depends only on dimension, it is clear 
that fs,k = ?^e,fc£~"^ is a bounded quantity, with j{e, k) = Cnf'^'^'^. 

We make use of a lemma, proved in [2] as Lemma 5.3, which says that for any set 
S confined to a ball of radius M with capacity 7, if if) is the capacitary potential of 
S and A^ the fundamental solution centered at 0, then 

Wx)--iN{x)\ < ^N(x) 
\x\ 

for all |a;| > ^. We reproduce the proof in an appendix for the reader's convenience. 
For our needs, we consider the behavior of this inequality at the "intermediate scale" 

represented by a^ = e^^^^, to be precise, for ^ag < \x\ < Aag, making the substitutions 

n— 1 / \ 1 

M -)■ Me^^, and 7 -^ 7(A;, e)e''-^: 

\ifj{x)-e''-^^{e,k)N{x)\ < Ce^^^+^^ 
We assume e is sufficiently small that Me "-2 < ^. 

13 



Consider now a single ball B^, with a set S^ representing the obstacle set, contained 
inside B u-i . Let ipei^) be the capacitary potential of S^ with respect to M", and 
let tfe(x) be the standard corrector as defined in §2, using the equivalent radius r^. 
Let ^/^(x) be the usual bridging function, which is smooth, takes the value of 1 inside 
Ba^, outside B2as, with iVr/^l < ^ and \Ari^\ < ^, and define our new corrector by 

We = 'ne{x)l/je{x) + (1 - r]e{x))We(x) 

repeating the construction abou each Xe,k- 

4.2 Corrector conditions 

We need to check that our new corrector satisfies the corrector conditions. That ([T]) 
and ([2]) are satisfied is easy, as is the L^ portion of ([3]): 

1. That we have We = 1 on Ss,k is clear from the definition of ips- 

2. That t&e — )■ in L'^{D) is also clear, since iWg — )■ in L'^{D), and We only 
differs from 10^ inside [J B2a^{xk), which, summed over all balls, has volume 

0(6 £ "-2 ) = 0{e 2(n-2))^ an(j ^s both functions are bounded by 1, this 
means the L^ norm of their difference goes to 0. 

To show that || VtDe||L2(£)) is uniformly bounded, and that (jl]) are satisfied, is a 
slightly more delicate matter. 
Let us begin by writing 

VWs = r/eV^e + (1 - Ve)^We + {ipe " We)Vr]e (7) 

We can certainly write 

|Vu;,|' < C (?7'|V^,|' + (1 - vf\Ww,\'' + (^, - «;.)'|V77,|2) 

The first two terms present no difficulty for us, as ifje has energy of order 0(e"'~^), 
which summed over 0{e^~^) balls is 0(1). Similarly, we already know that the energy 
of We is uniformly bounded by the arguments of the previous section. Thus the 
problem becomes one of estimating [■^^ — ly^l on the support of V?]^, or in the annulus 
i?2a. \ Ba^ , which wc settle with the following estimate. 

Lemma 9. Inside the annulus B4^a^ \ Bi^^, we have that 

li^e-Wel =0{e) 

Proof. Not surprisingly, we will be comparing We with ■y{e)e"'~^N{x), and 'j{e)e"'^^N{x) 
with ipe- We claim first that on the annulus, 

\we{x) - ^{£)£''-^N{x)\ = 0{e) 
14 



This is done by noting that 











n-2 

' £ 






7(^)5"- 


-^N{x) = 




\X "-2 


, and recalhng that 














We{x) 


^2-n 


— £^~" 






r2-n 
' £ 


— £^^" 


A bit of calculation 


gives 


US that 






1 




= e^-vr^ 


1 


1 \ 


^2—n p2~n 


v' 1- 


efr^) 



£"-Vr 0(5) = 0{f) 



The constant portion of u'e(x) is of order e can be easily checked. Now taking advan- 
tage of the fact that for a^ < \x\ < 2as, we have \x\ = 0(6 "^^), we see that 

\we - 7(£)£"-^iV(x)| = 0(5")^ + 0{e) = 0{e) 

The next part, to compare 'j{e)e'^~'^N{x) with '^/'^(x), uses our lemma, which tells us 
that 

|7(£)£"'iiV(x) - 4(x)| < C£^+^l^ = o{e) 

on the scale with which we are concerned. Thus, 



\We -i^el = 0{e) 



With this estimate in hand, it is not hard to see that 

\Vv£\'\ws-^£\'<C-^. 



u 



n(n-^) 



Integrating over the annulus, which has volume 0{a^) = 0{e "-2 ), we find that 



/ 



' B2a^\Ba^ 

Summing over the entire volume of 0(£^^") balls, we see that the contribution of this 
term goes to as e — )■ 0. 

An estimate related to Lemma [9] we will use later is the following: 

Lemma 10. Inside the annulus B2a^ \ Ba^, we have that 

Proof. We notice that inside the annulus B^a^ \ Bi^ , both ^^ and We are harmonic 
functions, and hence we can apply interior gradient estimates, the distance to the 
boundary being of 0{ae). D 
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4.3 Convergence of the integral term 

We are left only with showing that w satisfies the requirements of (jl]). We do this 
by arguing that the change in J^ ^Vw^ ■ Vv^dx introduced by replacing w^ by We is 
negligible. Integrating by parts on a single ball, we find that 



(fNWe ■ VVedx = / (pVed^WedSx — / VeV(p ' Vl&e + Vs(j)/!S.Wedx 

BeHD Jd{Bend) JB.nD 

We will argue that, as e — ?■ 0, both interior terms go to zero, and the boundary term 
splits into three parts: a purely positive component lying along S^ which is when 
Vs = on S , a. term lying on dD which tends to 0, and a term that recovers the limit 
we are looking for. 

The first interior term can be estimated directly: since i&e ^ in H^{D), and v^ 
is bounded in L?{D), it is clear that that term goes to as e — )■ 0. 

The second interior term relies on an estimate of Aty^. Inside Ba^ and B^r \ B2a^, 
We is equal to ipe and w^, respectively, and these are both harmonic functions. Thus, 
we need only consider the behavior of the interior term in the annulus i?2ae \ Ba^. 
Noting that 

Aw, = r]eAtlje + (1 - Ve)AWe + A?7e(^e - We) + 2Vr]e ■ {Vt/Je ' VtW,) 



we proceed to bound each term. Af/^, = Aws = 0, and by construction we have 

„ 2(n-3/2) ^ n-3/2 

lAr/el < ^ = 0{e "-2 ), \Vrje\ < ^ = 0{e "-2 ). We invoke Lemmas [9] and [10] to 
control the terms in ipe — We, the result being that 



Thus 



\Awe\ < C— = Ce^^ 
ai 



^VeAwe < Cale^^ = ©(e"-^/^) 



B.nD 



Repeating this integral over all balls, we find that the second interior term has con- 
tribution 0{e2), which goes to as e —t- 0. 

The boundary oi B^Cl D has two components: 

d{B, n D) = {dD n B,) U {dBe n D). 

On dBe n D, we have w, = We, and so we can apply the result from §3.2.21 to show 
that it recovers our desired limit. Thus, we are only concerned with showing that we 
can neglect the contribution of the other portion of the boundary, specifically, that 
portion of it lying in dD fl B^ . 
We seek to control 

(pVed^WedSr, 
dDnB^ 
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which is more profitably written as 



We consider each term individually. 

The arguments of §3.2.1l apply to /g^p,^ (1 — i]e)duU!edSx, and suffice to show that 
term is 0. 

For the term ^qj^^^ ("^e — Ws)duVedSx, we apply Lemma to control ^^ — w^ on 
the annulus where Vr^^ is supported, as well as the estimate \Vri^\ < — , to get 

-. n-3/2 

Since Vrj^ only lives on the annulus B2a^\Ba^, we can use the bound that \dDr\B2a^\ < 
Cd^~^ which gives us the basic estimate that 



{ipe - We)dur]edS^ 

dDriBe 



< Ce'+^-'/^ = o(e 



n~l\ 



Summing over 0{e^ ") such balls, we see that the contribution of this term goes to 
0. 

For the last term, we have to do a little more work: 

Lemma 11. 

(f}Ve'nedu'ipe,kdSx -^ 



?I 



f, JdDr\B^(xk) 

as e ^ 0. 

Proof. As usual, we do our proof over a single ball, showing that the term is so small 
that, when summed over all balls, it still goes to 0. We consider the nature of i/^^. 
It is well known that (see, for example, the classical paper of Littman, Stampacchia, 
and Weinberger [7] , sections 5 and 6) the capacitary potential ^ of a set 5* solves the 
equation 

AV^ = 0- 

where u is a measure supported on the boundary of S', known as the capacitary 
distribution, and (j{S) is the capacity of S. Furthermore, if) can be thought of as the 
Newtonian potential of a. For our purposes, if we let o"e be the capacitary distribution 
corresponding to the portion of the obstacle set living in the center of the ball S^, 
this means that 

Js, \x-yr 
Let us analyze the normal derivative of ._\n-2 - Since x and y both lie on the C^'" 
surface dD, \iy{x) — z/(?/)| < C\x — |/|°. We recall that 

^ 1 / ^ X — y 

Vi r^ = 2 - n). V 

\x — yy^ ^ \x — yy- 
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Thus, (2 — n) I ^_7^n ■ ^{x) is the term to be estimated. To estimate the normal deriva- 
tive, we set local coordinates centered at x, such that the tangent plane to dD at 
X is {xn = 0}, and dD is locally represented by a C^'" function x„ = r(a;'), where 
x' represents the first n — 1 components. Thus, in local coordinates, we have that 
\xn — Vnl < C\x — y|^^°. Interestingly, u{x) = e^ in our coordinates, so 



|x — 2/1" ^ Ix — y|" \x — yl" ^ " 

Returning to our original integral, we find that, since r]^ is only supported on a ball 
of radius 2ae, we have 

JdDnBe 



< / \(l)Vsriedu1pe\dS^ 

JdDnB2a, 

JdonBon, 



< Uv,\\l 

'dDnB2a. 



dDCiBia^ 



JdDnBon., J F ^ y| 



'dDnB2. 

From here, we use Fubini to change the order of integration, and integrate in polar 
coordinates w.r.t x, remembering that a is supported inside a radius much smaller 
than Og 

<CUv,\\lo^ fa'^daiy) 

<C||0t;,|Uoo7,£"-i+-'^ 

Since v^ is uniformly in L°°, and 7^ is bounded, we have that this term is o{e"'^^), 
which yields the desired result. D 
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